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Abstract. The upper bound on the value of the surface gravity, (ja, for neutron stars with equations of state 
respecting iisound < c, is derived. This bound is inversely proportional to the maximum allowable mass Mmax, and 
it reads Qb < 1.411 x 10^^ (Mmax/Afo)"^ cm s~'^. It implies an absolute upper bound 7.4 x 10^* cm s~^ if one uses 
the 2a lower bound on the neutron mass measured recently in 4U1700-37, I.QMq. A correlation between and 
the compactness parameter 2GM/Rc? for baryonic stars is analyzed. The properties of Qs of strange quark stars 
and its upper bounds are discussed using the scaling properties of the strange-star models. 
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1. Introduction 

Gravitational acceleration on the stellar surface, usually 
called the surface gravity, is an important parameter of 
the theory of stellar atmospheres, in particular the at- 
mospheres of neutron stars (see, e.g., Zavlin & Pavlov 
2002). For neutron stars, gs should be defined taking 
into account the space-time curvature. The value of gs ^ 
few X 10^'* cm s^^ is involved the relation between internal 
and surface temperature of neutron stars (Potckhin et al. 
2003 and references therein). 

The surface gravity of neutron stars is by many orders 
of magnitude larger than for other stars; it is ~ 10^ times 
larger than for white dwarfs, and 10^ times stronger than 
the gravity at the solar surface. As we show in the present 
paper, at a given neutron-star mass the value of gs de- 
pends very strongly on the largely unknown equation of 
state (EOS) of dense matter at supra-nuclear densities. It 
is therefore of interest to derive upper bounds on gg. In 
Sect. 121 we derive such upper bounds on gs resulting from 
the condition of subluminality of the EOS (speed of sound 
not exceeding c). We show how an upper bound on gs can 
be obtained from measured neutron star masses. In Sect.|31 
we analyze the EOS-dependence of gs for a set of 31 models 
of dense matter and compare the maximal surface gravity, 
reached at the maximum allowable mass, with subluminal 
upper bounds. In Sect. ^ we consider surface gravity of 
strange quark stars. We derive a scaling formula for gs 
of strange stars and use it to relate the maximum sur- 



face gravity for various models of such stars. Concluding 
remarks are presented in Sect. O 

2. Bounds on surface gravity 

It is well known that (see, e.g., Shapiro & Teukolsky 1983) 

(1) 



9s 



where M is the gravitational mass of the star and R 
is the stellar circumferential radius, and effects of rota- 
tion are neglected. For a "standard neutron star" with 
M = 1.4 Mq and i? = 10 km one has gs = 2.43 x 
10^** cm s~^. It is therefore convenient to measure gs in 
units of 10^^ cm s^^, and to use gs,ii = .9s/(10^^ cm s~^). 

Let us introduce the dimensionless compactness pa- 
rameter X = Tg/R = 2GM/Rc^ where rg is the 
Schwarzschild radius. Then 



5s,i4(a^,M) = 15.21 



X M 



(2) 
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Generally, physical theories have to be causal. In the 
literature, the condition of causality is usually replaced 
by the constraint of subluminality of the EOS, Usound = 
(dP/dp)i/2 < c (see, e.g., Rhoades & Ruffini 1974, Hartle 
1978, Kalogera & Baym 1996, Keiser & Polyzou 1996, 
Glendenning 1997, Koranda et al. 1997, Haensel et al. 
1999), but strictly speaking, subluminality is not equiv- 
alent to causality. For instance, Bludman & Ruderman 
(1968) constructed causal and Lorenz- invariant EOSs with 
^^sound > c. In practice, the difference seems to be unim- 
portant; the counterexamples of Bludman & Ruderman 
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(1968) correspond actually to a highly excited medium 
(Bludman & Ruderman 1970). Nevertheless, for the sake 
of rigor, we will call an EOS respecting the condition 
f^sound ^ c "subluminal" , and an EOS violating this con- 
dition will be called "superluminal" . 

For subluminal EOSs of dense matter, a strict upper 
bound on x is Xmax = 0.708 (see, e.g., Haensel et al. 1999 
and references therein). Using this value we get, from Eq. 
((21), an upper bound for a surface gravity of non- rotating 
neutron star of mass M, 



^ SL 
.9s, 14 ^ 5max.l4 



14.11 



M 



(3) 



The maximal mass Mmax obtained from the realistic 
EOS should be higher than the largest measured neutron 
star mass, M^f^. Known rotation frequencies of neutron 
stars are sufficiently small compared to the mass-shedding 
(Keplerian) limit, therefore we can use the non- rotating 
approximation for observed neutron stars. We thus have, 
in terms of gs, an upper bound based on the observations 
and subluminality conditions 



5s,i4 < 14.11 



(4) 



The higher M^^^, the stronger the constraints on the 
maximum allowable gs (the smaller the value of the upper 
bound). A precisely measured mass of the Hulse- Taylor 
pulsar, (1.4408 ± 0.0006) Mq (Weisberg & Taylor 2003), 
gives the upper bound 9.79 x 10^^ cm s^^. According to 
Clark et al. (2002), neutron star in the high-mass X- 
ray binary 4U1700-37 has, at the 2a confidence level, 
-^^obs > 1.9 Mq. This implies an upper bound lower than 
7.4 X lO^'' cm s~^. It should be mentioned, however, that 
the result of Clark et al. (2002) has to used with care, 
because of large errors. Moreover, because the compact 
object in 4U1700-37 is neither an X-ray pulsar nor an X- 
ray burster, it could in principle be a low-mass black hole. 

The expression on the right-hand-side of Q deserves 
an additional comment. Apart from the numerical con- 
stant, it is identical to the upper bound on the frequency 
of stable rigid rotation of neutron stars with subluminal 
EOS (Koranda et al. 1997). An alternative derivation of 
an approximate (but very precise) upper bound on rota- 
tion frequency, relating it to the maximum stellar com- 
pactness (surface redshift) for subluminal EOS, was given 
by Haensel et al. (1999). In the latter work, the starting 
point was the "empirical formula" for the maximum ro- 
tation frequency. Its upper bound was also obtained by 
maximizing a function of a; = fg/R. 

3. Equation of state of dense matter and surface 
gravity 

We have calculated gs for neutron-star models based on 
a set of 31 baryonic EOSs of dense matter. The EOSs 
were obtained under different assumptions regarding the 
composition of the matter at p ^ 2po, where po = 2.7 x 
10^^ g cm"'^ is the normal nuclear density. 
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Fig. 1. Surface gravity in the units of 10^ cm s~^ against 
neutron-star mass for 31 EOSs of dense baryonic mat- 
ter. Only stable configurations are shown, so that the 
curves terminate at the maximum allowable mass. Thick 
dashed line: maximally-stiff-core EOS. Remaining EOSs: 
thin lines. For further explanation see the text. 



The EOSs can be subdivided into several groups. 
Within one group of models, the matter consists of nucle- 
ons and leptons (Baldo et al. 1997, Bombaci 1995, Balberg 
& Gal 1997, Balberg et al. 1999, model I of Bethe & 
Johnson 1974 (usually called BJI), Pandharipande 1971, 
Pandharipande & Ravenhall 1989, Douchin & Haensel 
2001, Wiringa et al. 1988, Walecka 1974, Haensel et al. 
1980, Akmal et al. 1998). Within the second group, the 
matter is assumed to consist of nucleons, hyperons and 
leptons (Glendenning 1985, Balberg & Gal 1997, Weber 
et al. 1991). The third group involves an exotic high- 
density phase: de-confined quark matter mixed with bary- 
onic matter (Glendenning 1997), pion condensate (Muto 
& Tatsumi 1990) and kaon condensate (Kubis 2001). 
Finally, one EOS is the so called "maximally-stiff-core" 
(MSG) EOS. It consists of the BJI EOS below the baryon 
density 0.3 fm~^, matched continuously to the EOS with 
Wsound = c at higher density. 

In Fig. n we show the surface gravity versus stellar 
(gravitational) mass for all 31 baryonic EOSs. The values 
of gs for a given mass strongly depend on the EOS. For 
M — 1.44 Mq, 5s_i4 ranges from 1.1 to 2.8, while for 
M = 1.9 Mq the predicted values are between 1.4 and 
4.7. 

We considered also three EOSs of self-bound abso- 
lutely stable strange quark matter, forming hypothetical 
strange stars (Zdunik et al. 2000, Dcy et al. 1998) - they 
will be studied separately in Sect.^ 
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Fig. 2. Plots of (7s, 14 versus compactness parameter 
2GM/Rc^. Notation as in Fig. ^ Thick solid line rep- 
resents approximate formula, Eq. |(SJ). 
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Fig. 3. 
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Maximum values of surface gravity in units of 
(centers of the symbols at the bottom end of 
the vertical segments) for 34 baryonic and quark EOSs of 
dense matter, and the upper-bounds at the maximum al- 
lowable mass for these EOSs obtained from Eq. Q (dashes 
at the upper end of the vertical segments). Nucleonic EOS: 
B1,B2 - Baldo et al. (1997); D - Bombaci (1995); El, E3 
Balberg & Gal (1997), Balberg et al. (1999); C - model 
I of Bethe & Johnson (1974); J - Pandharipande (1971); 
F - Pandharipande & Ravenhall (1989); H - Douchin & 
Haensel (2001); LI, L2, L3 - Wiringa et al. (1988); - 
Walecka (1974); I - Haensel et al. (1980); A - Akmal et 
al. (1998). Hyperonic EOSs: G1-G5 - Glendenning (1985); 
E2, E4 - Balberg & Gal (1997), Balberg et al. (1999); Ml, 
M2 - Weber et al. (1991). EOSs with exotic high-density 
phase: G6-G8 - Glendenning (1997); K - Kubis (2001); 
PI, P2 - Muto & Tatsumi (1990). Strange quark matter: 
SI, S2 - Zdunik et al. (2000); S3 - Dey et al. (1998). N 
labels the MSG (maximally-stiff-core) EOS. For further 
explanations see the text. 



Gontrary to the significant scatter of the gs{M) plots 
for baryonic stars, gs{x) is much less EOS-dependent: it 
can be rather well (within better than 20%) reproduced 
by the approximate formula 



9sM 



5a;5/4 



(5) 



Actually, the precision of this formula is much higher if 
we exclude the superluminal EOSs and the unrealistically 
stiff ones (too high incompressibility of nuclear matter at 
saturation) which yield Mmax ^2.5 M0. Putting then the 
upper bound Xmax — 0.708 into ((SJ yields an approximate 
"realistic upper bound" of 6 x 10^^ cm s^^ on gs of neutron 
stars. 

The maximum value of gs for stable stars with a 
given EOS is reached at the maximum allowable mass. 
The values of (/s,max for the selected EOSs are shown in 
Fig. 131 For comparison, we show also the absolute up- 
per bounds, Eq. |(2Jl, at Mmax- The dense matter mod- 
els were divided into groups denoted by specific sym- 
bols. Subluminal models involving only nucleons (o) give 
5s,max ~ (3 - 5) X 10" cm s-2. The stiffer the EOS, the 
closer gs,max to the subluminal upper bound f/^'^^xl-^max)- 

Subluminal hyperonic EOSs (•) give s's.max which is 
typically lower than 4 x 10^^ cm s~^. For these EOSs, 
fl's.max can be as small as one-fifth of the upper bound 

5max(-^rnax)- 

The subluminal EOSs with an exotic high-density 
phase (©) have relatively low 5s. max- A phase transi- 
tion softens the EOS, lowering the radius. Simultaneously, 
however, the softening leads to a decrease of Mmax- The 
latter effect dominates over the former one. This strongly 
pushes up the upper bound gs(-^max)- If the stellar inte- 
rior consists mostly of a mixed quark-baryon phase (EOSs 



2 X 10 cm s , only 



tn ri^^ 
ys,max 



G6, G7, G8 in Fig.EJ, then ^s 

one fifth of 5^L^(Mmax)- 

The MSG EOS (o) yields 5s, max which is quite close 
A similar situation occurs for the EOSs that 
give Minax models with superluminal cores, labeled 
Their values of 5s, max range from 4 x 10"'^'* cm s~^ to 
6 X lO^** cm s~^. Note that even for those EOSs 5s, max < 
5s max- This is because the compactness parameter is al- 
ways smaller than the upper bound for subluminal EOSs, 
a;(Mi„ax) < a;max = 0.708, as is shown in Fig.|2| 

4. Surface gravity of strange stars 

The gs{M) dependence for strange stars, presented in the 
upper panel of Fig. ^ is very different from that of bary- 
onic stars, due to a different R{M) dependence. The model 
dependence of gs{M) for strange stars is very strong: at 
1.44 Mq, the value of 5s, 14 ranges from 2.0 to 5.5. However, 
this strong EOS dependence can be explained in terms of 
the scaling properties of strange star models. 

EOSs of self-bound quark matter (strange matter) are 
derived from various different models of the quark struc- 
ture of hadrons. Despite differences between the underly- 
ing models, the EOSs relevant for stable models of strange 
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Fig. 4. Surface gravity in units of 10^^ cm versus 
strange star mass (upper panel) and compactness parame- 
ter (lower panel), for three EOSs of strange quark matter. 
Dashed line in the lower panel is obtained by transforming 
(see text for details) the S2 EOS curve into S3 EOS curve 
using Eq. ij?)). 



stars can be well represented (fitted) by a linear relation 
between the pressure P and the mass density p (Zdunik 
2000), 



P = ac^{p- ps) ■ 



(6) 



The parameters of stellar configuration, calculated for an 
EOS given by Eq. ©, are connected to those obtained for 
an EOS with p'^ ^ ps by the scaling relations (Witten 1984, 



Haensel et al. 1986, Zdunik 2000). For example, points on 
the M[R) curves are related by 



M = I ^ 



M' 



R' 



(7) 



Therefore, at a fixed a, the ratio M/R (and therefore x = 
Tg/R) does not depend on ps, and the maximum surface 
gravity of strange stars scales as 



— 9s 



(8) 



The values of a for the EOSs SI, S2, and S3 range 
within 0.30 < a < 0.46. We have the exact ratio 
5s.max(S3)/(/s,max(S2) — 1.81, whilc the scaling factor 
[ps(S3)/ps(S2)]i/2 = 1.63. The large difference in maxi- 
mum surface gravities can therefore be accounted for by 
scaling with respect to ps', in this context, the a depen- 
dence is sufficiently weak to be of neglected in the sim- 
plest approach. As we see in the lower panel in Fig. 21 the 
precision of this scaling increases with decreasing strange 
star compactness. The scaling becomes very precise for 
X < 0.2. The explanation of this behaviour is simple: the 
density within such strange stars is nearly constant, so 
that the dependence on a is negligible. 

5. Discussion and conclusion 

We derived an upper bound on the surface gravity of 
static neutron stars with subluminal EOSs. Maximum gs 
is reached at the maximum allowable mass Afmax- Even at 
the shortest observed pulsar period of 1.56 ms the effect 
of rotation on the M^ax is very small, so that the static 
approximation is justified. The upper bound on is in- 
versely proportional to the highest measured neutron-star 
mass M^f'' and is 7.4 x lO^"' cm s~^ if M^f^ is replaced by 
1.9 M0, which is the 95% confidence-level lower bound on 
mass of the compact star in 4U1 700-37 measured by Clark 
et al. (2002). However, this upper bound on gs should be 
taken as highly unreliable. Firstly, Clark et al. (2002) used 
a specific model of the companion star and of the binary to 
evaluate the compact-star mass. Secondly, one cannot ex- 
clude, unfortunately, that the compact object in 4U1700- 
37 is actually a low-mass black hole and not a neutron 
star (Clark et al. 2002 give arguments for and against the 
black- hole presence). 

We have studied gs for a set of 31 EOSs of baryonic 
matter. The dependence of on the stellar mass M is very 
sensitive to the EOS. On the contrary, the dependence of 
gs on the stellar compactness rg/ R has a generic character 
for baryonic EOSs. The maximum surface gravity ^7s,max 
is sensitive to the EOS of dense matter and ranges from 
about 2 X 10^"' cm s"^ to 5 x 10^** cm s for subluminal 
baryonic EOSs. 

The dependence of on the mass and compactness of 
strange stars is very different from that of baryonic stars. 
However, the range of gs,max = (3 — 6) x 10^^ cm s~^ is 
quite similar to baryonic stars. 
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We hope that the EOS-sensitive features of gs will be 
useful in extracting information about the EOS of dense 
matter, for instance, by combining the values of gs ob- 
tained fitting the thermal component of the observed pho- 
ton spectra with atmospheric models and the surface red- 
shift measured for the identified spectral lines. 

Acknowledgements. This work was motivated by a question of 
A. Majczyna concerning the possible range of surface grav- 
ity in neutron star atmospheres, which is gratefully acknowl- 
edged. We are grateful to D.G. Yakovlev for the reading of 
the manuscript and for helpful remarks. This work was par- 
tially supported by the KBN grants no. 2P03D.019.24 and 
5P03D.020.20. 

References 

Akmal, A., Pandharipande, V.R., & Ravenhall, D.G. 1998, 

Phys. Rev. C, 58, 1804 
Balberg, S., Gal, A. 1997, Nucl. Phys. A, 625, 435 
Balberg, S., Lichtenstadt, I., & Cook, G.B. 1999, ApJS, 121, 

515 

Baldo, M., Bombaci, I., & Burgio, G.F. 1997, A&A, 328, 274 
Bethe, H.A., & Johnson, M.B. 1974, Nucl. Phys. A, 230, 1 
Bludman, S.A., & Ruderman, M.A., 1968, Phys. Rev., 170, 
1176 

Bludman, S.A., & Ruderman, M.A., 1970, Phys. Rev. D, 1, 
3243 

Bombaci, I. 1995, in Perspectives on Theoretical Nuclear 
Physics, ed. I. Bombaci, A. Bonaccorso, A. Fabrocini et 
al., 223 

Clark, J. S., Goodwin, S. P., Crowther, P. A., Kaper, L., 
Fairbairn, M., Langer, N., Brocksopp, C, 2002, A&A, 392, 
909 

Dey, M., Bombaci, I., Dey, J., Ray, S., & Samanta, B.C., 1998, 

Phys. Lett. B, 438, 123 
Douchin, F., & Haensel, P. 2001, A&A, 380, 151 
Glendenning N. K. 1985, ApJ, 293, 470 

Glendenning, N.K. 1997, Compact Stars: Nuclear Physics, 
Particle Physics and General Relativity (Springer, New 
York) 

Haensel, P., Kutschera, M., & Proszynski, M. 1980, A&A, 102, 
299 

Haensel, P., Lasota, J. P., & Zdunik, J. L., 1999, A&A, 344, 
151 

Haensel P., Zdunik J. L., & Schaefer R., 1986, A&A, 160, 121 
Hartle, J.B., 1978, Phys. Repts., 46, 201 
Kalogera, V., Baym, G., 1996, ApJ, 470, L61 
Keiser, B., Polyzou, W., 1996, Phys. Rev. C, 54, 2023 
Koranda, S., Stergioulas, N., Friedman, J.L., 1997, ApJ, 488, 
799 

Kubis, S. 2001, Ph.D. Thesis, Jagiellonian University (unpub- 
lished) 

Pandharipande, V.R. 1971, Nucl. Phys. A, 147, 641 
Pandharipande, V.R., & Ravenhall, D.G. 1989, in Proc. NATO 

Advanced Research Workshop on nuclear matter and heavy 

ion collisions, Les Houches, 1989, ed. M. Soyeur et al. 

(Plenum, New York, 1989), 103 
Potekhin, A.Y., Yakovlev, D.G., & Chabrier, G., 2003, ApJ, 

594, 404 

Rhoades, C.E., & Ruffini, R., 1974, Phys. Rev. Lett., 32, 324 



Shapiro, S. L., Teukolsky, S. A., 1983, Black holes, white 
dwarfs, and neutron stars: The physics of compact objects. 
New York, Wiley-Interscience 

Muto T., Tatsumi T., 1990, Prog. Theoret. Phys. 83, 499 

Walecka, J. D., 1974, Ann. Phys., 83, 491 

Weber, F., Glendenning, N.K., & Weigel, M. K., 1991, ApJ, 
373, 579 

Weisberg, J.M., Taylor, J.H., 2003, in: Radio Pulsars, edited 
by M. Bailes, Nice D.J., & Thorsett S.E., ASP Conference 
Series (to be published); e-print astro-ph/0211217 

Wiringa, R.B., Fiks, V., & Fabrocini, A. 1988, Phys. Rev. C, 
38, 1010 

Witten, E., 1984, Phys. Rev. D, 30, 272 

Zavlin, V.E., & Pavlov, G.G., 2002, in: Proceedings of the 
270. WE-Heraeus Seminar on Neutron Stars, Pulsars, 
and Supernova Remnants. MPE Report 278. Edited 
by W. Becker, H. Lesch, and J. Truemper. Garching 
bei Muenchen: Max-Planck-Institut fuer extraterrestrische 
Physik, p. 263 

Zdunik, J.L., 2000, A&A, 359, 311 



